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ABSTRACT 

We present our photometric study of the stellar association NGC 602 in the wing of the Small Magellanic 
Cloud (SMC). The data were taken in the filters F555W^ and .F814W using the Advanced Camera for Surveys 
(ACS) on-board the Hubble Space Telescope (HST). Photometry was performed using the ACS module of 
the stellar photometry package DOLPHOT. We detected more than 5,500 stars with a magnitude range of 
14 m555 ~ 28 mag. Three prominent stellar concentrations are identified with star counts in the observed 
field, the association NGC 602 itself, and two clusters, one of them not being currently in any known catalog. 
The Color-Magnitude Diagrams (CMDs) of both clusters show features typical for young open clusters, while 
that of the association reveals bright main sequence (MS) and faint pre-main sequence (PMS) stars as the 
members of the system. We construct the initial mass spectrum (IMS) of the association by applying an age- 
independent method of counting the PMS stars within evolutionary tracks, while for the bright MS stars we 
transform their magnitudes to masses with the use of mass-luminosity relations. The IMS of NGC 602 is 
found to be well represented by a single-power law, corresponding to an Initial Mass Function (IMF) of slope 
r w — 1.2 for 1 ^ M/M 45. This indicates that the shape of the IMF of a star forming system in the SMC 
for stars with masses higher than 1 M Q seems to be quite similar to the field IMF in the solar neighborhood. 
Subject headings: galaxies: star clusters — Magellanic Clouds — open clusters: individual(NGC 602) — stars: 
evolution — stars: pre-main-sequence stars: luminosity function, mass function 



1. INTRODUCTION 

The Small Magellanic Cloud (SMC), the second closest 
undisrupted neighboring dwarf galaxy to our own (after the 
Large Magellanic Cloud; LMC) is an ideal laboratory to in- 
vestigate star formation at an environment similar to the early 
universe. Its metallicity is by a factor of 5 lower than the solar 
metallicity (Luck et al. 1998; Ven 1999), and its dust-to-gas 
ratio is measured to be up to 10 times lower than in the Milky 
Way (e.g. Stanimirovic et al. 2000). Therefore, the SMC may 
be considered as an excellent local template for studying pri- 
mordial star formation, providing insight into the processes at 
work in the early universe. The recent star formation process 
in the SMC is characterized by a rich sample of H II regions 
(Henize 1956; Davies et al. 1976) linked to young stellar sys- 
tems, the stellar associations, known for their early-type stel- 
lar content and loose structure (e.g. Kontizas et al. 1999). The 
small distance of the SMC from us (« 60 kpc) allows us to re- 
solve from space individual stars down to My « 10 mag, and 
consequently the study of its low-mass young stellar content 
and the corresponding Initial Mass Function (IMF) becomes 
unusually accessible. 

Stellar associations contain the richest sample of young 
bright stars in a galaxy. Our knowledge on the young mas- 
sive stars of the Magellanic Clouds (MCs) has been collected 

1 Research supported by the Deutsche Forschungsgemeinschaft (German 
Research Foundation) 



from ground-based studies of young stellar associations (see 
e.g. Massey 2006 and references therein). However, our 
knowledge of the low-mass stellar membership and the cor- 
responding IMF in star-forming regions of the MCs is quite 
incomplete. A first attempt to define the low-mass population 
of a stellar association in the MCs was made by Gouliermis 
et al. (2005). In their photometric study of the LMC associa- 
tion LH 52 with data from the Wide-Field Planetary Camera 
2 (WFPC2), these authors test the hypothesis that sub-solar- 
mass stars can be detected in the MCs. They found that all 
faint main sequence stars, which are observed in the area of 
the association, belong to the general field of the LMC and not 
to the system. The field-subtracted mass function of LH 52, 
which accounts for its IMF was found by these authors with a 
slope T ~ —1.1 for main sequence stars down to M ~ 1 M©. 
This slope is comparable, but somewhat more shallow than a 
typical Salpeter (1955; T 1.35) IMF. An IMF well repro- 
duced by a power law with a slope consistent with Salpeter' s 
was also found by Sirianni et al. (2002) in the central regions 
(within 30") of the young SMC star cluster NGC 330. 

The investigation by Gouliermis et al. (2005) revealed, for 
the first time, that the low-mass stellar content of stellar as- 
sociations in the MCs is still in its pre-main sequence (PMS) 
phase, as is the case of nearby galactic associations (Briceno 
et al. 2007). Gouliermis et al. (2006a) discovered in LH 52 
about 500 low-mass PMS stars down to ~ 0.5 M Q , easily 
distinguishable in the 771555 — m 555 Color-Magnitude 
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FIG. 1. — Color-composite image of the field around NGC 602, constructed from observations with HST/ACS in the filters F555W, F814W and 7^658^. 
The association itself lies at the center of the field, located almost in the middle of the shell-like emission nebula N 90. Directly north of the association 
lies the open cluster B 164 (Bruck 1976), and the newly discovered "cluster A" covers the north-western edge of the field. All three systems are confirmed 
as statistically significant stellar concentrations from star counts shown in Figure [4] Credit: NASA, ESA, and the Hubble Heritage Team (STScI/AURA)- 
ESA/Hubble Collaboration. 



TABLE 1 
Log of the observations 



Visit 


RA 


DEC 


Tcxpo 




Dataset 


(J2000.0) 


F555W 


F814W 


J92F05 


01 29 28.21 


-73 33 16.7 


430 


453 




01 29 28.52 


-73 33 16.4 


430 


453 




01 29 27.89 


-73 33 16.9 


430 


453 




01 29 27.57 


-73 33 17.1 


430 


453 




01 29 27.57 


-73 33 17.1 


3 


2 




01 29 27.25 


-73 33 17.4 


430 


453 


J92FA6 


01 29 27.57 


-73 33 17.1 


3 


2 


Note. 


— Datasets refer to HST archive catalog. 


Exposure times (T cxpo 



per filter are given in seconds. Units of right ascension are hours, minutes, 
and seconds, and units of declination are degrees, arcminutes, and arcsec- 
onds. 



Diagram (CMD), extending the sub-solar stellar membership 
of MCs stellar associations to their PMS populations. Natu- 
rally, these stars should be considered as the best tracers of the 
sub-solar IMF in the MCs, but the WFPC2 data, being lim- 
ited by completeness, allowed the construction of the PMS 
IMF of LH 52 only for the mass range 0.8 - 1.4 M . This 
IMF was found to correspond to a power-law with a slope 
of r ~ —1.26. Previous PMS studies in the MCs focus on 
the surrounding field of supernova 1987A (e.g. Panagia et al. 
2000), where the sample is limited, and the star-burst of 30 
Doradus (e.g. Zinnecker 1998; Sirianni et al. 2000; Brand- 
ner et al. 2001), where crowding and high extinction limit the 
detection of stars to 1 - 2 M Q (Romaniello et al. 2006). 

Over the last years, stellar associations in the SMC gained 
increasing scientific interest, as the investigation of their 
bright OB stars was complimented by new studies on their 
PMS stars (e.g. Nota et al. 2006; Gouliermis et al. 2006b; 
Sabbi et al. 2007; Hennekemper et al. 2008). These studies 
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TABLE 2 

Sample from the photometric catalog of stars found in this 
study in the region of ngc 602/n 90 with hst/acs imaging 





RA. 


DECL. 


# 


(J2000.0) 


(J2000.0) 


1 


01 29 54.828 


-73 32 31.524 


2 


01 29 24.581 


-73 33 16.236 


3 


01 29 14.220 


-73 31 53.256 


4 


01 29 33.310 


-73 33 44.496 


5 


01 29 31.051 


-73 33 43.056 


6 


01 29 28.301 


-73 31 57.432 


7 


01 29 36.485 


-73 34 12.108 


8 


01 29 32.837 


-73 33 39.024 


9 


01 29 31.992 


-73 33 30.924 


10 


01 29 27.768 


-73 32 59.496 


11 


01 29 18.286 


-73 32 33.468 


12 


01 29 31.094 


-73 33 43.272 


13 


01 29 31.411 


-73 33 42.336 


14 


01 29 14.249 


-73 31 53.220 


15 


01 29 32.942 


-73 34 10.956 



m 555 


CT 555 


?71gl4 


^"814 


(mag) 


(mag) 


(mag) 


(mag) 


14.690 


0.004 


13.217 


0.003 


13.881 


0.002 


14.128 


0.003 


14.717 


0.003 


13.354 


0.002 


14.270 


0.002 


14.508 


0.004 


15.371 


0.004 


15.612 


0.007 


15.341 


0.004 


15.288 


0.006 


16.412 


0.007 


14.993 


0.005 


15.467 


0.004 


15.673 


0.007 


15.594 


0.005 


15.796 


0.007 


15.730 


0.005 


15.699 


0.007 


16.052 


0.006 


15.322 


0.006 


15.816 


0.005 


16.074 


0.008 


15.663 


0.005 


15.871 


0.008 


17.910 


0.020 


16.655 


0.018 


16.232 


0.006 


16.377 


0.010 



NOTE. — Magnitudes are given in the Vega system. Units of right ascen- 
sion are hours, minutes, and seconds, and units of declination are degrees, ar- 
cminutes, and arcseconds. Tablets available in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 



TABLE 3 

Populations Summary in the Identified Stellar Systems 



Name ristars p (pc 2 ) / PM s /ums /lms 

NGC 602 855 5.6 0.86 0.07 0.04 

B 164 356 2.9 0.03 0.27 0.62 

Cluster A 171 2.0 0.05 0.21 0.63 

Field 230 LI 0.16 0.08 0.66 

NOTE. — The total stellar numbers of the systems, n s tars, are field- 
subtracted (see § 4.3.1), and refer to the area of each system defined with 
star counts. The total number of field stars, shown in the last row, refers to 
the area defined as the best representative of the field due to its emptiness. 
The fractions / give the number of one stellar species over the total stellar 
numbers in every area. The most significant numbers are highlighted. It is 
interesting to note that both B 164 and Cluster A have almost identical num- 
bers of UMS and LMS stars, and that their LMS fractions are comparable to 
that of the field. 

were facilitated by the unique combination of high-resolving 
power with a wide field-of-view provided by the Advanced 
Camera for Surveys (ACS) on-board HST. In the present 
study we deal with the IMF of the young stellar association 
NGC 602, related to the bright H II region LHA 115-N90, 
or in short N 90 (Henize 1956). The region NGC 602/N 90, 
being located in the wing of the SMC, has the advantage of 
avoiding the densest parts of the galaxy. As a consequence, 
the confusion of its stellar content with that of the general 
SMC field, which extends in a wide depth along the line-of- 
sight, is rather low. McCumber et al. (2005) found in their 
investigation on field stars in the SMC wing that after a long 
period of constant star formation, there was a burst in activ- 
ity within the last 1 Gyr. The lack of 7-ray, X-ray or far-UV 
sources in the region of NGC 602/N 90 suggests that this ex- 
cess in the star formation of the area was possibly triggered 
by encounters with the LMC and/or the Milky Way. Most 
likely, all major star-forming events in the SMC are initiated 
by these tidal interactions. In the first part of our investiga- 
tion of the region NGC 602/N 90 (Gouliermis et al. 2007) 
we present the results of our photometric study from obser- 
vations with the Infrared Array Camera (IRAQ on-board the 
Spitzer Space Telescope. We report the detection of 22 candi- 



date Young Stellar Objects (YSOs), which are located at the 
edge of N 90, suggesting that they are the products of star 
formation triggered by the central association NGC 602. 

The recent release of the HST/ACS observations of the re- 
gion of NGC 602/N 90 offers a unique opportunity for the 
comprehensive photometric analysis of both its bright and 
faint stellar content, and the construction of the complete IMF 
throughout the whole observed mass range. Our photome- 
try, which reaches the limit of TO555 <; 26.5 mag, providing 
one of the most complete stellar samples ever collected for 
star-forming regions in the SMC, is described in § 2. We in- 
vestigate the morphology of the observed region and the spa- 
tial distribution of the detected stars in § 3, and in § 4 we 
discuss the nature of the detected stellar species and present 
the CMDs of the three major stellar concentrations in the ob- 
served field. In § 4 we also present the subtraction of the 
contribution of the field from the observed CMDs, and we 
estimate the reddening towards the identified systems. A de- 
tailed description of the stellar members of the association 
NGC 602 is also given in § 4. The Luminosity Functions of 
the identified stellar systems are constructed in § 5. In § 6 a 
comprehensive construction of the IMF of NGC 602, as well 
as of the two additional star clusters in the region is presented. 
Conclusions are given in § 7. 

2. OBSERVATIONS AND PHOTOMETRY 

The data used in our analysis are obtained with the Wide 
Field Channel (WFC) of the ACS within the HST GO Program 
10248 on July 14 (Dataset: J92F05) and July 18 (Dataset: 
J92FA6) 2004. Images are taken in the F555W and F8UW 
filters. A single ACS/WFC pointing covering w 3f4x3'4 was 
centered on NGC 602/N 90. At the distance of the SMC this 
field-of-view (FoV) corresponds to an extent of « 57x57 pc. 
The first visit (J92F05) includes five long exposures taken 
with a dithering pattern in order to cover the inter-chip gap 
of the camera and one short exposure with a single pointing 
to cope with the saturation of the brightest stars. Within the 
second visit (J92FA6) a short exposure was also taken in or- 
der to remove cosmic rays that were unintentionally included 
during the first visit. In TableQ]we give a summary of the dif- 
ferent data sets used in this investigation. A color-composite 
image of the observed field is shown in Fig.[T] 

Photometry was performed with the ACS module of the 
package DOLPHOT 2 (Ver. 1 .0), especially designed for ACS. 
We followed the photometric process, as it is described by 
Gouliermis et al. (2006b). The pipeline-reduced FITS files 
were obtained from the HST Data Archive. We used the pack- 
age multidrizzle (Koekemoer et al. 2002) following the ACS 
data handbook to clean the images of residual warm pixels 
and cosmic rays and to construct a deep reference drizzled 
image to be used for the photometry. All exposures were pho- 
tometered simultaneously, using the ^8141^ drizzled frame 
as the position reference. Photometric calibrations and trans- 
formations were made according to Sirianni et al. (2005). 
Charge Transfer Efficiency corrections were applied accord- 
ing to ACS ISR 04-06. We cleaned our photometric catalog 
from bad detections using DOLPHOT's star quality parame- 
ters, as described also by Gouliermis et al. (2006b). The final 
photometric catalog includes 5,626 stars in total detected in 
both filters. The full length of this catalog with the RA, DEC 

2 The ACS module of DOLPHOT is an adaptation of the photom- 
etry package HSTphot (Dolphin 2000). It can be downloaded from 

http : //pur cell . as . arizona . edu/dolphot / 
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FIG. 2. — Typical photometric uncertainties in both ms55 and mgi^ bands 
as derived by DOLPHOT from all data sets. 
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FIG. 3. — Completeness for all stars in the observed FoV and in the area 
selected as the most representative of the local field the SMC. The field be- 
ing less contaminated by the nebula and less crowded than the association 
NGC 602 exhibits a higher completeness than the overall FoV. 

coordinates (in J2000) and their magnitudes in each filter (in 
the Vega system) of these stars is available in electronic form 
in Table [2] Astrometric solutions for the detected stars are 
derived from the F814W^ drizzled frame with the use of the 
application xy2sky of the WCSTools package 3 . 

Typical uncertainties of our photometry as a function of 
the magnitude for both filters are given in Fig. |2] The com- 
pleteness of the data was evaluated by artificial star experi- 
ments with the use of lists of almost 400,000 artificial stars 
created with the utility acsfakelist of DOLPHOT. The com- 
pleteness is found to be spatially variable, depending on the 
crowding of each region. This is shown in Fig. [3] where the 
completeness factors for the whole observed field, as well as 
of the area selected as the most representative of the back- 
ground field, due to its emptiness, are plotted for both fil- 
ters. The area, which from here on we refer to as simply 
the "field", is selected as the most empty area of the west- 
ern part of the observed FoV. This selection was based on our 
recent results from SpitzerflRAC, which show that this part 
of NGC 602/N 90 is the less contaminated by dust emission 
(Gouliermis et al. 2007). 

3 T he WCSTools package is available at C fA at 
|http : //tdc-www .harvard, edu/soft ware/west ools/| 
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FIG. 4. — Isodensity contour map for all stars with photometric errors 
ay < 0.1 mag using grid elements of size 7", shown as the small box at 
the bottom-left corner of the map. The first isopleth corresponds to the mean 
background density. The next isopleths are plotted in steps of the standard 
deviation a. The bluish contours correspond to stellar density equal or higher 
than the 3cr level. All three stellar systems in the region are revealed as sig- 
nificant stellar concentrations. 

3. MORPHOLOGY OF THE REGION OF NGC 602/N 90 

The prominent bright stellar concentration almost at the 
center of the observed ACS/WFC field, shown in Fig. [TJ is 
the association NGC 602. The bright emission nebula, N 90, 
can also be seen in the image, because of the Ha emission, 
included in the F658N filter (centered at 658.4 nm). A small 
contribution from [Nil] is also included in this filter. The as- 
sociation is surrounded by a ring-like feature highlighted by 
the gas emission, and ribbon-like structures pointing towards 
the center can be seen to the southeastern and northwestern of 
the ring. These structures have been identified by Gouliermis 
et al. (2007) as the SMC analogs of the "Pillars of creation" 
originally detected in the Galactic H II region M 16 (Hester et 
al. 1996). They clearly indicate ongoing star formation on the 
rim of the ring-like structure, probably driven by the stellar 
winds and UV radiation of the brightest stars in NGC 602. 

The stellar association is not the only stellar concentration 
inside the observed FoV. North of the association lies the open 
cluster B 164 (Bruck 1976) which is also known as NGC 602b 
(Westerlund 1964). Next to B 164 and to the west lies a third 
stellar concentration. It shows to be a cluster, which does not 
belong to any catalog of known objects in this region. From 
here on we refer to this cluster as "Cluster A". All three stellar 
systems are easily revealed by the spatial distribution of all 
stars detected with our photometry in the region, as discussed 
in the next section. 

3.1. Spatial Distribution of the Stars in NGC 602/N 90 

In order to reveal the spatial distribution of the detected 
stars and the morphology of any stellar concentrations in the 
observed region, we performed star counts on our photometric 
catalog (a sample of this catalog is shown in Table [2]). Stars 
with photometric uncertainties ay > 0.1 mag were not con- 
sidered. The constructed isodensity contour map is shown in 
Fig-El The coordinate system used for this map, as well as the 
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ones shown later is the Cartesian (X,Y) system of the pixel co- 
ordinates of WFC, transformed to seconds of arc, in respect 
to the central pixel of the FoV. 

Star counts were performed in a quadrilateral grid divided 
in elements with sizes 140 x 140 WFC pix 2 each, which cor- 
respond to 7"x7" (or w 2 pc x 2 pc at the distance of the 
SMC). This size was selected as the optimum to map any fine 
structure, which might exist within the clusters. The selection 
of larger size for the grid elements would reveal a smoother 
version of the map of Fig. [4] The first isopleth in the con- 
tour map indicates the mean background density, whereas the 
subsequent levels are drawn in steps of la, where a is the 
standard deviation of the background density. All contours 
with density equal or higher than the 3cr level, which is cho- 
sen as the density threshold of the statistically significant stel- 
lar concentrations are drawn with bluish colors. All three star 
clusters in the region are revealed in this map with NGC 602 
being the dominant stellar system, seen almost at its center. 
The density peaks for both B 164 and Cluster A do exceed the 
3<t limit as well. Smaller compact density peaks can be seen 
in the map below the 3rr limit. They are typical density fluc- 
tuations, due to the small size selected for the grid elements. 

4. STELLAR POPULATIONS IN THE REGION OF NGC 602/N 90 
4.1. Color-Magnitude Diagram 

The TO555 — 771814, 7TJ555 color-magnitude diagram (CMD) 
of all stars in our photometric catalog is shown in 
Fig. [5] which indicates that our photometry provides ac- 
curate magnitudes (with 17555 $ 0.1 mag) for stars with 
771555 26.5 mag. Different groups of stellar types can eas- 
ily be distinguished in this CMD. There is a prominent well 
populated main sequence (MS), which extends from the de- 
tection limit up to TO555 » 14 mag, and a rich population of 
pre-main sequence (PMS) stars located on the red part of the 
MS. We select regions of the CMD, each covering different 
types of stars. The tentative limits of these regions are also 
shown in Fig. [5] We split the MS into its upper (brighter) and 
lower (fainter) parts (UMS and LMS respectively). Accord- 
ing to the ZAMS isochrone from the model grid designed by 
Girardi et al. (2002) for the ACS filter system, the limit be- 
tween UMS and LMS stars of 771555 = 22.6 mag corresponds 
to a stellar mass of M w 1.2 M Q . 

The brighter limit of the PMS region in the CMD is be- 
ing chosen around 771555 « 22.3 mag. Probably, a brighter 
limit would allow us to include the brighter PMS stars and 
the transition region between PMS and MS (turn-on), but this 
part of the CMD is contaminated by the old stellar population 
of the field of the galaxy (e.g. McCumber et al. 2005), as 
seen by its loose red giant branch (RGB). The low numbers 
of widely scattered stars in the RGB part of the CMD and the 
absence of any clearly detectable turn-off indicates that only a 
small fraction of the evolved field stars can be seen in the ob- 
served region, probably due to extinction by the dust. Indeed, 
Spitzer/IRAC observations of NGC 602/N 90 have shown in 
the 8 /im band strong dust emission from the area, which sur- 
rounds the central "hole" in the H II region, where NGC 602 
is located (Gouliermis et al. 2007). 

4.2. Stellar Populations in the Systems of the Region 

Isodensity contour maps from counts of the three most im- 
portant stellar types in the CMD of Fig. [5] exhibit a strong 
spatial distinction between UMS, LMS and PMS stars. This 
is shown in the bottom panels of Fig. [6] (Figs. [6^, 03, |6}i and 
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FIG. 5. — The 1TI555 — mgi4, 711555 CMD of all stars detected in the region 
of NGC 602/N 90. Stars with photometric uncertainties larger than 0. 1 mag 
(Fig. [2} are shown as grey dots. The drawn lines indicate the tentative limits, 
which separate the different stellar species included in the catalog of detected 
stars: The Upper Main Sequence (UMS), the Lower Main Sequence (LMS), 
and the Pre-Main Sequence (PMS) stars. 

|6^). The contour maps of this figure are constructed with the 
same method as the map of Fig. [4] (see § 3.1), and the low- 
est density isopleth shown corresponds to 3cr above the local 
background density. These maps reveal again all three impor- 
tant stellar systems in the region (NGC 602, B 164 and Clus- 
ter A), but they demonstrate that these systems are dominated 
by different types of stars. 

All three systems are detected by star counts of UMS stars. 
Although B 164 and Cluster A are found to contain signif- 
icant samples of LMS stars, they are not revealed at all in 
the map of the PMS stars (Fig. |6ji). On the other hand, the 
faint stellar populations of the association NGC 602 are cer- 
tainly PMS and not LMS stars. This demonstrates that both 
B 164 and Cluster A seem to be evolved small clusters, but 
NGC 602 is a younger stellar association, with a significant 
number of PMS stellar members. Consequently, in order to 
reveal the extent of each system based on its dominant stellar 
populations we construct their contour maps from star counts 
of these populations alone, shown in the top panels of Fig. [6] 
The map constructed for B 164 and Cluster A from the UMS 
and LMS stars is shown in Fig. [(J, while the one for NGC 602 
(from counts of the UMS and PMS stars) is shown on Fig. [6f. 
We define the boundaries of each stellar system from the 3a 
isopleth in these maps. The size of each system is defined by 
a rectangle, which contains this isopleth. 

We show in Table[3]the numbers of stars counted within the 
measured limits of each system (after the field contribution 
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FIG. 6. — Isodensity contour maps for the region of NGC 602/N 90 from star counts of different stellar populations (LMS, UMS, PMS; see §4.1), as they are 
revealed in the observed CMD. The maps of the whole region (bottom panels) have a grid resolution of 12", while the enlarged ones (top panels) of 5". In all 
maps the first drawn density isopleth corresponds to the 3<r significance. 



subtracted, see 4.3.1), the surface density and the correspond- 
ing fractions of PMS, UMS, and LMS stars included in these 
limits. These fractions show a clear concentration of PMS 
stars in NGC 602, while Cluster A and B 164 clearly host 
the largest amount of LMS stars. Although, both Cluster A 
and B 164 show to host larger numbers of UMS stars, it is 
NGC 602 that hosts the brightest MS stars in the region, as it 
will be shown later from the CMDs of the systems. Consid- 
ering that the limits for the UMS stars in the CMD have been 
selected to include stars down to 771555 ~ 22.5 mag, the larger 
fraction of UMS stars in Cluster A and B 164 is due to their 
fainter UMS members. 

4.3. The Star Clusters B 164 and Cluster A 

According to our selection of different stellar types in the 
observed CMD (§ 4.1), about 55% of the stars in our catalog 
are LMS stars. Fig. [6^ shows that these stars are mostly con- 
centrated in B 164 and Cluster A, but in a rather loose manner 
(there are only three isopleths above the local background). 
On the other hand, the concentration of UMS stars in these 
clusters shows a peak at 9a above the local background in 
B 164 (Fig. [6J5). The detailed contour map of the area around 
these clusters, shown in Fig. [6};, shows the prominent con- 
centration of both UMS and LMS stars in both of them (and 
mostly in B 164) indicating their evolved nature. The absence 
of PMS stars in both clusters as seen in the map of Fig. [6}i, 
certainly affects their individual CMDs. In order to construct 
these CMDs we consider the contribution of the stellar pop- 
ulation of the general field of SMC in the region, as it is ob- 
served within the observed ACS pointing. We apply a Monte 
Carlo method for the subtraction of this contribution from the 
CMD observed in the specific area of each cluster. 

4.3. 1 . Subtraction of Contaminating Field Stars 

The typical procedure for the field subtraction of a CMD 
with the use of the Monte Carlo method divides the CMD in a 
grid and performs a comparison between the numbers of stars 



found in each grid element of the observed CMD of the clus- 
ter and the one of the region selected as representative of the 
general field. Although this method is efficient for compact 
populous star clusters, it might produce unwanted artifacts for 
smaller not compact clusters or associations. Therefore, our 
Monte Carlo method for the field subtraction of the observed 
CMDs of the clusters slightly differs from the typical proce- 
dure. Our field-subtraction procedure does not use a prede- 
fined grid on the whole extent of the CMDs, but an individual 
subregion of the CMD is considered around every single star 
observed in the area of the cluster. Each such subregion has 
the shape of an ellipse with semi-axes A(V — I) = 0.15 mag 
and AV = 0.5 mag. 

The corresponding numbers of field and cluster stars inside 
each of these ellipses (np for the field and «c f° r the cluster) 
are counted, and corrected for incompleteness according to 
the completeness factors estimated per magnitude range for 
each of the considered areas. The number of field stars per 
ellipse expected to be present in the CMD of the cluster, rip, 
is given by the number of field stars, np, normalized to the 
surface of the area of the cluster. The probability that the 
central star of each ellipse is actually a field star is given by 
n F /ric- If e.g. np = then all the stars inside the ellipse 
are cluster stars and the probability is 0, while if n F /nc > 1 
then we can assume that this star is actually a field star. For 
intermediate values, < ri F /nc < 1, it is randomly decided 
whether to account this star for a field star or not. 

We applied this method to decontaminate the CMD of the 
area of each system from the contribution of the field. The 
field-subtracted CMDs of both Cluster A and B 164 are shown 
in Fig. [7] Their most remarkable feature is the clear MS, 
which extends from the detection limit up to the brightest stars 
in the clusters, with magnitudes 771555 ~ 16.8 mag for Clus- 
ter A and 777,555 ~ 15-5 mag for B 164. An older limit for 
the age of both clusters can be estimated based on their bright 
MS stars, with isochrone fitting. We used the models by Gi- 
rardi et al. (2002), especially designed for the filter system of 
ACS, and we found an age of « 160 Myr for Cluster A, while 
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Cluster A. The peaks of the distributions were estimated assuming that they 
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E(B-V) = 0.06±0.02 mag, for Cluster AE(B-V) = 0.04±0.02 mag, 
and for B 164 E(B - V) = 0.07 ± 0.02 mag. These values show that all 
three clusters suffer from rather low optical extinction. 



B 164 is found be not older than w 80 Myr. It should be noted 
that if the only bright red star seen in the CMD of Cluster A 
is also considered, an age of r ~ 80 Myr may also represent 
this cluster. If indeed the two clusters have similar ages, and 
considering that they are very close to each other, one may 
speculate that both clusters form a double system, product of 
the same star formation event. However, the available data are 
not sufficient to support or dismiss this hypothesis. 

In Fig. [7] both CMDs are shown with the corresponding 
isochrones overplotted. In these CMDs, as well as in all 
CMDs shown here we assume a distance modulus (j,q ~ 
18.9 ± 0.1 mag (Dolphin et al. 2001). The largest observed 
stellar masses for MS stars in each cluster can be estimated 
from the evolutionary models considered. These masses are 
- 4.3 M Q for Cluster A and - 5.7 M for B 164. 



4.3.2. Interstellar Reddening 



For the isochrones plotted in Fig. [7] no specific interstel- 
lar reddening was considered. The comparison of the loci of 
the individual UMS stars to the models allowed us an accu- 
rate determination of the reddening within the area of each 
cluster. The corresponding distributions of the stars in Clus- 
ter A and B 164 according to their reddening are shown in 
Fig. [8] (the distribution for NGC 602 is also shown for com- 
parison; see § 4.4). These plots show that a mean color excess 
of E(B - V) = 0.04 ± 0.02 mag can be considered for Clus- 
ter A and E(B - V) — 0.07 ± 0.02 mag for B 164. These 
values, as well as the one for NGC 602, are in excellent agree- 
ment with previous reddening measurements, which vary be- 
tween E{B — V) ~ 0.03 and 0.07 mag for clusters located 
in the Wing of the SMC (Piatti et al. 2007), and between 
E(B — V) ~ 0.05 and 0.09 mag for southern fields of the 
galaxy (Noel et al. 2007). They are also in good agreement 
with previous measurements on SMC intermediate-age clus- 
ters (e.g. Alcaino et al. 2003; Hunter et al. 2003; Rochau 
et al. 2007). The reddening values found for Cluster A 
and B 164 correspond to an extinction of Ay w 0.1 mag 
and Ay w 0.2 mag respectively assuming an extinction law 
A v = R V E(B - V) with R v = 3.1 (e.g. Koornneef 1983). 

4.4. 77ie Stellar Association NGC 602 

As mentioned earlier (§ 4.2) NGC 602 does not host sig- 
nificant numbers of LMS stars (Fig. |6^), but rather a large 
concentration of PMS stars (Fig. [(J). The system is also char- 
acterized by a prominent UMS stellar population (Fig. [6}i). 
Fig. |6f shows in more detail the morphology of the system 
based on the spatial distribution of both UMS and PMS stars. 
This map shows clear signatures of sub-clustering within the 
system (as is the case for Cluster A and B 164, shown in 
Fig. [6}:), which is demonstrated by the clumpy behavior of 
the density in the map. Although, this phenomenon could be 
due to the selection of a finer grid resolution in the star counts 
(with a grid element size of ps 5'.'0 instead of 12'.'0, which was 
used for the larger maps), from the positions of the UMS stars 
we verified that they are indeed concentrated in smaller sub- 
regions within NGC 602, the most prominent being the one 
located to the south-east (bottom-right) of the map in Fig . [6f. 
This subgroup has also an over-density of PMS stars located 
with an offset of w 5" in respect to the density peak of the 
UMS stars. A second important subgroup of UMS stars is 
located w 10" east of the first, and as shown in Fig. |6f it is 
somewhat elongated. 

The CMD of the region of NGC 602 is shown in Fig. in 
comparison to the CMD of the region of the field. NGC 602 
shows a more well populated and brighter UMS than Clus- 
ter A and B 164, clearly suggesting its youthfulness. We iden- 
tified the 10 brightest stars of NGC 602 as early-type with 
spectral types between 05.5 and B0. 5. These "photometric 
spectral types" are consistent (±0.5 spectral type) with those 
previously derived by Hutchings et al. (1991) from optical 
and UV spectroscopy. Moreover, the most prominent feature 
in its CMD is the large number of stars in the PMS part, which 
is completely empty in the CMDs of both Cluster A and B 164 
and of the field 4 . For the subtraction of the contribution of 

4 The selection of the region of the field could not be made without includ- 
ing some PMS stars. The reason is that no control field was observed away 
from NGC 602 and therefore we selected the most "empty" and remote part 
of the observed FoV as the best representative of the field. The inclusion of 
PMS stars in the field may lead to the underestimation of the the number of 
PMS stars in NGC 602. However, as derived from the numbers of Table 3, 
this cannot affect more than 5% of the total PMS members of NGC 602. 
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FIG. 9. — The 777555 — 777814, "1555 CMD of the stars in the area of the as- 
sociation NGC 602 (without any subtraction of the contribution of the field), 
and of the area selected as the most representative of the general background 
field of the SMC. In the latter, features which are characteristic of the SMC 
wing stellar population, such as a turnoff at rj 22 mag (McCumber et al. 
2005), corresponding to the overplotted isochrone of age ~ 1 Gyr, can be 
distinguished. 

the field population in the CMD of NGC 602 we applied the 
same Monte Carlo method (described in § 4.3.1) as for Clus- 
ter A and B 164. The derived "clean" CMD of the stars of 
NGC 602 alone, shown in Fig. [To] clearly indicates that all 
the observed LMS stars belong to the field and not the associ- 
ation. The distribution of interstellar reddening for the UMS 
stars of NGC 602 (Pig. [HJ also shows very low extinction. The 
mean color excess obtained assuming a Gaussian distribution 
is E(B -V) = 0.06 ± 0.02 mag. 

In the CMD of Fig. [TO] it can be seen that the positions of 
the PMS stars show a spread as if they belong to groups of 
different ages. The observed broadening of the loci of the 
PMS stars in the CMD implies a wide span in ages of over 10 
Myr. Three indicative isochrones of ages 2, 4 and 10 Myr are 
overplotted in Fig.[l0]to demonstrate this spread. If the CMD 
broadening of the PMS stars is produced by an age-spread 
or not is very important in understanding the star formation 
history of the region. The existence of an age-spread would 
suggest that these stars did not originate from a single star 
formation event, but they were formed over a longer period of 
time. On the other hand if all PMS stars are the product of a 
single star formation event, then there should be other factors 
that affect their CMD positions. 

The importance of the age of PMS stars lies on the fact 
that these stars are the live records of star formation that took 
place in a star forming region within the last ^ 30 Myr (e.g. 
Bricefio et al. 2007). Consequently, if we assume that differ- 
ent age-groups of PMS stars can be easily distinguished by 
their locations in the CMD, they should be the best tracers 
of any sequential star formation event that might have taken 
place within their host stellar systems. The age-spread of 
PMS stars in young stellar systems of our Galaxy has been 
discussed by e.g. Palla (2005), who notes that the star for- 
mation history of the ONC extends long in the past, although 
at a reduced rate, and by Preibisch & Zinnecker (2007), who 
found no evidence for an age dispersion in Upper Scorpius 
OB association, although small age spreads of w 12 Myr can- 
not be excluded. 

However, considering that PMS stars have a rather peculiar 
behavior, as there are several effects that change their position 
in the CMD, the retrieval of any signature of age-spread in a 
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FIG. 10.— The 771555 - m 8 i4, m 55 5 CMD of the stars in NGC 602, 
after the contribution of the field has been statistically subtracted with the 
Monte Carlo method. A clear upper MS (down to 771555 — 21.5 mag), as 
well as a sequence of faint PMS stars (with 777,555 ~ 21.5 mag) can be easily 
distinguished to represent the stellar populations of the association. The loci 
of the PMS stars show a prominent broadening. This makes the estimation 
of the age of these stars rather difficult. A discussion on the origins of such 
a spread in the positions of PMS stars in the CMD is presented in § 14.41 
Characteristic isochrones from the grid of PMS models by Siess et al. (2000) 
are overplotted. The birth-line and the ZAMS from the same grid of models 
are also plotted. 

PMS population with single-epoch photometry alone is not a 
trivial task. Hennekemper et al. (2008) using the same ob- 
servational material of another SMC association, NGC 346, 
showed that there are various physical properties of the PMS 
stars of the system, which affect their loci in the CMD. This 
produces a broadening, which could be misinterpreted as an 
age-spread. These properties, typical for T Tauri stars of our 
Galaxy, are unresolved binarity, variability and differential 
extinction. Although it is most certain that these characteris- 
tics are responsible for most of the observed CMD broadening 
of the PMS stars (Hennekemper et al. 2008), one cannot ex- 
clude the possibility that there might be also a true age-spread 
among these stars, hidden in the CMD of NGC 602. We cur- 
rently explore this possibility in another paper by modeling 
the observed CMD with the construction of synthetic CMDs 
for NGC 602 (M. Schmalzl et al., in preparation). 

5. LUMINOSITY FUNCTIONS 

We constructed the luminosity function (LF) of all three 
stellar systems by counting the stars in magnitude bins of 
0.5 mag in both 771555 an d m-814- The m555-LFs of the systems 
are shown in Fig. [TTJfor the field subtracted stellar samples 
shown in Figs. [7]and[T0] The LFs of the counted numbers of 
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FIG. 1 1 . — The r?i555 -luminosity functions (LFs) for all three stellar sys- 
tems. The limit of the 50% completeness is indicated in each LF by the dotted 
lines. The measured numbers of stars are shown with the dashed lines, while 
the thick lines show the LFs after the stellar numbers have been corrected 
for incompleteness. In all shown LFs the contribution of field stars has been 
statistically subtracted. 

stars are shown with dashed lines, while the LFs shown with 
thick lines are corrected for incompleteness. The dotted ver- 
tical line for every plot marks the 50% completeness limit for 
each system. The LFs of both Cluster A and B 164 are very 
similar, and they both suffer from lower number statistics than 
the LF of NGC 602 shown in the bottom panel of Fig. QT] 
In the latter one can identify a small dip at 771555 ~ 22 mag. 
This feature has been theoretically predicted as a direct conse- 
quence of the stellar evolution of young star clusters, since the 
stars at 771555 ~ 21 mag gain brightness very fast, depleting 
thus this part of the CMD (Stahler & Fletcher 1991; Fletcher 
& Stahler 1994). 

6. MASS FUNCTION 

6.1. Mass Spectrum 

The stellar Initial Mass Function (IMF), £(logM), is the 
distribution of the stars in a stellar system according to their 
masses, by the time of their formation. The IMF gives the 
number of stars of the system per unit logarithmic (base 10) 
mass, dlogM, per unit area. There are several functional 
forms proposed to represent the IMF (see Kroupa 2002 for 
a review), but its intermediate- and high-mass part (down to 
M ~ 1 Mq) is generally characterized by a single-power law 
of the form £(logM) ~ M . Typical logarithmic derivative 
is r k —1.35, as estimated for stars in the solar neighbor- 
hood with masses 0.4 < M/M < 10 (Salpeter 1955). If 
the stars are counted in unit mass intervals, dM, per unit area, 
then their Initial Mass Spectrum (IMS), f(M), is constructed. 
IMS is generally described also by a power law of the form 



f(M) ~ M a . The IMS and IMF are connected by the rela- 
tion (Miller & Scalo 1979): 

/(M) = C(logM)^ (1) 

Consequently, the relation for their indexes (slopes) is: 

a = T - 1 (2) 

For the construction of the IMF or IMS a mass luminos- 
ity (ML) relation based on isochrone models is required for 
the translation of the observed stellar luminosities to masses. 
However, as shown in § 4.4, the loci of the PMS stars in the 
CMD of NGC 602 are well spread, producing a broadening 
which covers a wide range of isochrones (Fig.fTOl). Since the 
observed CMD broadening of PMS stars is not solely due to 
age differences, but also due to the characteristics of these 
stars (Hennekemper et al. 2008), a construction of a unique 
ML relation is rather impossible. Specifically, the observed 
optical magnitudes of young PMS stars are affected by several 
factors, especially the extinction/emission produced by the ac- 
cretion disk and/or envelope (e.g. Herbst et al. 2002; Sherry 
et al. 2004; Briceno et al. 2005). At ages less than 10 7 yr 
all low-mass PMS stars are variable, and the strength of this 
photometric variability depends on the star's magnetic activ- 
ity, the amount of circumstellar matter, and the accretion rate. 
Classical T Tauri stars are actively accreting gas from circum- 
stellar disks, producing irregular photometric variability with 
amplitudes up to 3 mag in the V band (Herbst et al. 1994). 
Naturally, such variations affect the estimation of the lumi- 
nosities and the effective temperatures for these PMS stars, as 
well as of their individual circumstellar and differential red- 
dening. Therefore, any attempt to fit isochrone models to the 
observed CMD and to determine the ages to these stars will 
be biased 5 . 

As a consequence a ML relation cannot be used for the 
determination of the masses of the observed PMS stars in 
NGC 602. Therefore, in the following sections we study the 
IMS of the association, as we constructed it by applying an 
age-independent approach of counting the PMS stars within 
PMS evolutionary tracks (e.g. Gouliermis et al. 2006a). For 
this method we use two different grids of evolutionary tracks, 
by Palla & Stahler (1999) and Siess et al. (2000), for stellar 
masses between 0.1 and 6 M Q . Both grids cover an age range 
starting at the birth-line up to ps 20 Myr, and they are avail- 
able in the standard Johnson-Cousins photometric system. In 
order to compare our CMD to these models we transformed 
our ACS 777555 an d 777814 magnitudes to the standard V and / 
respectively with the use of the transformations provided by 
Sirianni et al. (2005; § 8.3, and Appendix D). In Fig. [12] we 
show the transformed V - I, V CMD of NGC 602 with the 
case of counting PMS stars between tracks from the models 
by Siess et al. (2000). Although the lowest metallicity consid- 
ered for this grid of models is Z = 0.01 (~ 2.5 times higher 
than that of the SMC) it fits our data better than the one from 
Palla & Stahler (1999), which assumes the metallicity of the 
SMC of Z = 0.004. Both grids agree for low mass stars up to 
« 2.5 M Q . For higher masses the models by Palla & Stahler 
(1999) show a shift to the blue. This can be demonstrated 
by the comparison of the ZAMS from both grids of models 
(Fig. [T2l . The most massive part of the IMS, for stars with 
Al <; 6 Mq, is occupied by evolved MS stars, and therefore 

5 See Hennekemper et al. (2008; § 3) for a discussion on the limitations in 
the age estimation for the PMS stars in the SMC association NGC 346. 
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another grid of models for such stars, was used. These models 
were compiled by Girardi et al. (2002) for the standard photo- 
metric filter system and the metallicity of the SMC. It should 
be noted that the ZAMS of Siess et al. (2000) coincides bet- 
ter with the one from Girardi et al. (2002), which match very 
well the observed MS of the system. 

For the construction of the IMS we count all stars with 
a completeness > 50%, and therefore only stars with 
m 55 5 24.5 mag, corresponding to M jSi 1.1 M Q are 
counted. Considering that each evolutionary track corre- 
sponds to a specific mass, we constructed the IMS by count- 
ing the PMS stars between evolutionary tracks based on their 
loci in the CMD. We constructed thus, a histogram with bins 
defined by the available PMS tracks. It should be noted that 
building a histogram is the simplest, but by no means opti- 
mum, way of analyzing the distribution of masses. Given the 
uncertainties in our photometry, produced by the character- 
istics of the PMS stars, which translate into corresponding 
mass uncertainties, it is possible that some of the objects in 
every bin might actually have masses slightly outside the bin 
range. Moreover, Mafz Apellaniz & Ubeda (2005) found sig- 
nificant numerical biases in the determination of the slope of 
power laws from uniformly binned data using linear regres- 
sion. These biases are caused by the correlation between the 
number of stars per bin and the assigned weights, and are es- 
pecially important when the number of stars per bin is small. 
This result implies the existence of systematic errors in the 
values of IMFs calculated in this way. As an alternative, those 
authors proposed to use variable-size bins and divide the stars 
evenly among them. Such variable-size bins yield very small 
biases that are only weakly dependent on the number of stars 
per bin. 

Taking into account the above, in order to avoid artificial 
features in the distribution of the masses, we smoothed our 
counts by folding them with an adaptive kernel function using 
a simple boxcar average procedure of 3 mass-bins width. This 
process smoothes out any random fluctuations due to the fact 
that some stars could belong to neighboring bins due to the 
photometric uncertainties and inaccuracies of the evolution- 
ary tracks. Furthermore, in order to properly eliminate any 
numerical biases in the estimation of the IMS slope due to 
the bin sizes, the stellar numbers per mass-bin, N, which are 
used for the IMS, are the counted incompleteness-corrected 
stellar numbers, N, but normalized to a bin-size of 1 M Q , by 
dividing them to the sizes AM of their corresponding bins. 

6.2. The Initial Mass Spectrum ofNGC 602 

The constructed IMS of NGC 602 with the use of evolu- 
tionary tracks from both Palla & Stahler (1999) and Siess et 
al. (2000) is shown in Fig. Q~3] Small changes in the slope 
of this IMS, which are apparent in the figure, especially the 
"knee"-like feature at « 2.5 M Q , imply that maybe a multi- 
power law is appropriate to represent this IMS. However our 
modeling of the IMS with a weighted linear fit proved that a 
single-power law of the form N ~ M a describes best this 
IMS. The slopes derived for the IMS constructed with the use 
of models by Palla & Stahler (1999) and Siess et al. (2000) 
show that the use of both grids gives more or less the same re- 
sults within the derived uncertainties. Specifically, the derived 
slope for the IMS constructed with the use of the tracks by 
Palla & Stahler (1999) is found to be a ~ -2.0 ±0.3, and the 
slope derived with the use of the tracks by Siess et al. (2000) 
a ~ -2.4 ± 0.2 for the mass range 1 M/M Q ^ 45. It 




0.0 0.5 1.0 1.5 2.0 
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FIG. 12. — Construction of the Initial Mass Spectrum (IMS) of NGC 602 
by counting stars between evolutionary tracks. The PMS tracks by Siess et 
al. (2000) are shown in this plot. The ZAMS of both Siess et al. (2000) and 
Palla & Stahler (1999) are shown to demonstrate their differences for stars 
with V < 20 mag. The ZAMS from the grids of models by Girardi et al. 
(2002), which is found to fit very well the one by Siess et al. (2000) is plotted 
for the upper MS of the CMD. This ZAMS has been used for the construction 
of a mass-luminosity relation for the bright MS stars with M ^ 6 Mq for 
the derivation of their masses. Different PMS tracks and the stars included 
between them are plotted in different colors to be easier distinguished. 

should be noted that if we compare the derived IMS slopes 
for specific narrower mass ranges, the largest difference be- 
tween the used models is seen in the lower mass range of 
1 Af/M < 2.5 (although these differences are covered by 
the fitting errors). This is, however, expected from the great 
differences between the ZAMS of Palla & Stahler (1999) and 
the one of Siess et al. (2000) shown in Fig. [TO] which affects 
mostly stars with M ^ 2.5 M Q . 

Taking into account the results from the application of the 
linear regression for the IMS constructed with the use of both 
libraries of evolutionary tracks, the IMS of NGC 602 is best 
described by the function 

n ~ ikr 2 - 2±0 ' 3 (3) 

for the whole observed mass range of ~ 1 M Q up to w 45 M Q . 
The derived mean slope has a remarkable coincidence with 
the slope of a = —2.35 found by Salpeter (1955) for the so- 
lar neighborhood, and with the results of Kroupa (2002) who 
suggests a slope of a = —2.3 ± 0.3 for the mass range of 
0.5 < MMq < 100 as the average field star IMF. The IMS 
of Kroupa (2002) is represented by the shaded region over- 
plotted at the top panel of Fig. [13] 

A third IMS to be considered for comparing the IMS of 
NGC 602 is by Scalo (1998) with a slope of a = -2.7 ± 
0.5 for 1 < M/M n < 10 and a = -2.3 ± 0.5 for 
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FIG. 13. — The Initial Mass Spectrum (IMS) of NGC 602 as constructed 
by counting PMS stars between evolutionary tracks and with the use of a ML 
relation from the models by Girardi et al. (2002) for the brightest MS stars. 
The derived IMS is shown in each panel from counting PMS stars between 
tracks from both the models of Siess et al. (2000) (top IMS) and Palla & 
Stahler (1999) (bottom IMS). The shades at the top panel correspond to a 
Kroupa (2002) IMS with a = -2.3 ± 0.3, while the ones at the bottom 
panel to a Scalo (1998) one with —2.7 ^5 a J$ — 2.3. Stellar numbers TV 
are corrected to a bin-size of 1 Mq. A slope of a ~ —2.2 ± 0.3 is found to 
represent very well the IMS as it is constructed with the use of all considered 
grids of models. 



10 < M/Mq < 100. The slopes we find for the IMS of 
NGC 602 are comparable also with this IMS for M > 10 M . 
For smaller masses the slopes are more different, but still con- 
sistent within the errors. The IMS of Scalo (1998) is repre- 
sented by the shaded region overplotted at the bottom panel 
of Fig. [13] It is worth noting that our single-power law slope 
for the IMS of NGC 602 agrees very well with the one found 
by Gouliermis et al. (2006a) for the stellar association LH 52 
in the Large Magellanic Cloud, a = —2.3 ± 1.0. These au- 
thors could estimate the IMS slope only for a narrow mass 
range (0.8 ^ M/M Q ^ 1.4) due to observational limita- 
tions. A multi -power law was found by Preibisch et al. (2002) 
for the IMS of the whole Upper Sco OB association with the 
use of data from 2dF, ROSAT, and Hipparcos. These authors 
found a slope of a = -2.8 ± 0.5 for 0.6 < M/M < 2 and 
a = -2.6±0.5for2 < M/M Q < 20. Both these slopes are 
systematically steeper than the one we found for NGC 602. 



6.3. Mass Spectra of Cluster A and B 164 

As it is seen earlier from their CMDs, both B 164 and Clus- 
ter A are open clusters without any significant PMS popu- 
lation. Therefore, the construction of their mass spectra is 
more straightforward, since it can be done with the use of 
mass-luminosity (ML) relations. These relations can be de- 
rived from certain models for evolved MS stars (Girardi et al. 
2002), and they provide an accurate transformation of lumi- 
nosities into masses for the stars of these clusters. Naturally, 
the choice of the appropriate isochrones, as the most represen- 
tative for the age of the systems, is crucial for this procedure. 
For this selection we used the upper age limits found earlier 
(§ to be w 160 Myr for Cluster A and » 80 Myr for 
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FIG. 14. — The mass spectrum of both B 164 and Cluster A, constructed 
using ML-relations based on the evolutionary models by Girardi et al. (2002). 
Points derived from incomplete data (with completeness < 50%) are shown 
with arrows. The mass spectrum of both clusters is well represented by a 
single-power law comparable to a Salpeter IMS. 

B 164. After obtaining the masses of the cluster members we 
constructed the corresponding IMS of each cluster, based on 
the results of Mafz Apellaniz & Ubeda (2005), by counting 
the stars into variable-size mass-bins. 

A single-power law is also found to represent the IMS 
of both clusters for the whole observed mass range with 
M ?^ 1.5 M Q . The best-fitting slopes of the corresponding 
IMS are found to be: 

N ~ M- 2 - 2±0 - 2 for B 164, N ~ M~ 2 - 9±0 ' 5 for Cluster A 

(4) 

The IMS slope of Cluster A is found significantly steeper than 
that of B 164. Specifically, it is interesting to note that the 
IMS of B 164 is more similar to the one by Kroupa (2002) 
(a = -2.3 ± 0.5), whereas the IMS of Cluster A fits better 
to the one by Preibisch et al. (2002) (between a ~ —2.8 and 
-2.6 ±0.5) or by Scalo (1998) (a = -2.7 ± 0.5). 

7. CONCLUSIONS 

We present our photometric study from data obtained in 
the filters F555W and F8UW with HST/ACS in wide-field 
mode of the star forming region NGC 602/N 90 in the wing of 
the SMC. Our aim is to take advantage of the high-resolution 
efficiency of ACS for the construction of the IMF of the as- 
sociation NGC 602 for the whole observed stellar population 
of the system. We use the ACS module of the photometric 
package DOLPHOT, which is especially designed for imag- 
ing with the ACS and we provide the full photometric cata- 
log of all stars detected with short and long exposures in one 
ACS/WFC field that covers a region of about 3.'4 x 3f4 cen- 
tered on the young stellar association NGC 602. The region 
is found to host a mixture of stellar populations. Star counts 
revealed that there are three prominent stellar concentrations 
in the observed field: The association NGC 602 itself at the 
center, and two young open clusters, B 164 (Bruck 1976), and 
one currently un-cataloged, which we name "Cluster A", both 
located at the northern edge of the observed area. The CMDs 
of all three systems are contaminated by the stellar popula- 
tion of the general field of SMC in the region. We selected 
the most empty (and less contaminated from dust emission, 
as seen with Spitzer; Gouliermis et al. 2007) part of the ob- 
served area as the most representative of the general SMC 
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field. We, then, applied a Monte Carlo method for decontam- 
inating the CMDs of all three identified stellar systems from 
the contribution of the field. 

The CMDs of B 164 and Cluster A show features typical for 
young open clusters, with a fully populated main sequence ex- 
tending from the bright part (with only a couple of bright stars 
in each cluster) down to the detection limit of 771555 ~ 28 mag. 
By comparing the brightest MS stars of both B 164 and Clus- 
ter A we were able to confine their maximum age to be around 
80 Myr and 160 Myr respectively. If we consider one red 
super-giant in the CMD of Cluster A as cluster member, then 
both clusters have the same age and their close separation sug- 
gests that they were possibly formed during the same star for- 
mation event. On the other hand, NGC 602 is a much younger 
system with some bright nebulosity still in its vicinity. The 
CMD of the association reveals that the stellar members of 
NGC 602 belong to a rich upper bright main sequence 6 down 
to TO555 ~ 21 mag and to a prominent red sequence of faint 
pre-main sequence stars down to the detection limit. In gen- 
eral, the observed sequence of PMS stars corresponds to an 
age of roughly 4 Myr. However, there is also a significant 
amount of PMS stars that can be fitted by a ~ 10 Myr model. 
The "field", the area selected as the most representative of 
the general SMC field, is characterized by a well populated 
low main sequence and a rather poor red giant branch with a 
turnoff at 777,555 ~ 22 mag. Few stars brighter than the turnoff, 
and few PMS stars were also included in this area, since it is 
located within the observed FoV. 

We constructed the luminosity functions of all three stel- 
lar systems, by counting the stars in magnitude bins of 0.5 
mag in both 771,555 and 777.814. Concerning the mass spectrum 
of the association, the PMS stars in the CMD of NGC 602 
show a prominent broadening, which does not allow a direct 
measurement of the masses of these stars through a single 
ML-relation. This broadening might be the result of an age 
spread, but it could also be due to the characteristics of the 
PMS stars, such as variability and extinction. In any case the 
accurate measurement of their age and the estimation of their 
masses becomes a rather difficult task. Therefore, we con- 
structed the initial mass spectrum of the association with an 
age-independent method based on counting the PMS stars be- 
tween evolutionary tracks. We used the grids of models by 
both Palla & Stahler (1999) and Siess et al. (2000). For the 

6 The 10 brightest stars in the system are identified as early-type with spec- 
tral types between 05.5 and B0.5 



bright main sequence stars we used a ML-relation obtained 
from the evolutionary models by Girardi et al. (2002). The 
IMS of the association is found to be well represented by a 
single-power law, corresponding to an IMF of slope T « — 1.2 
for 1 < M/Mq < 45, similar to the field IMF in the so- 
lar neighborhood (Salpeter 1955; Kroupa 2002). No signifi- 
cant differences are found between the IMF derived with the 
use of Palla & Stahler (1999) models and the one from the 
models by Siess et al. (2000), although these grids of models 
are found to be quite different from each other for stars with 
20 mag. We also constructed the mass spectra of 
the clusters B 164 and Cluster A, and we found that the mass 
spectrum of B 164 corresponds to a mass function similar to 
NGC 602 with V w -1.2, while the one of Cluster A is found 
to be steeper with T « — 1.9. 

The study of NGC 602/N 90 with observations from 
Spitzer/IRAC by Gouliermis et al. (2007) showed evidence 
of ongoing star formation on the rim of the H II region. Can- 
didate YSOs were identified by these authors, located around 
the region of NGC 602, and their formation is most likely trig- 
gered through a sequential event initiated at the central associ- 
ation. This hypothesis is currently under investigation within 
our study of the recent star formation history of this region 
through the reconstruction of the observed sequence of PMS 
stars in and around the association in synthetic CMDs (M. 
Schmalzl et al., in preparation). 
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